The Stochastic Dynamics of an Array of Atomic Force Microscopes in a Viscous Fluid 
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We consider the stochastic dynamics of an array of two closely spaced atomic force microscope 
cantilevers in a viscous fluid for use as a possible biomolecule sensor. The cantilevers are not 
driven externally, as is common in applications of atomic force microscopy, and we explore the 
stochastic cantilever dynamics due to the constant buffeting of fluid particles by Brownian motion. 
The stochastic dynamics of two adjacent cantilevers are correlated due to long range effects of the 
viscous fluid. Using a recently proposed thermodynamic approach the hydrodynamic correlations are 
quantified for precise experimental conditions through deterministic numerical simulations. Results 
are presented for an array of two readily available atomic force microscope cantilevers. It is shown 
that the force on a cantilever due to the fluid correlations with an adjacent cantilever is more than 
3 times smaller than the Brownian force on an individual cantilever. Our results indicate that 
measurements of the correlations in the displacement of an array of atomic force microscopes can 
detect piconewton forces with microsecond time resolution. 



I. INTRODUCTION 

The advent of micro and nanotechnology has ushered 
forth measurement techniques with unprecedented sen- 
sitivity 0, S Si 0- For example the atomic force mi- 
croscope (AFM) |5| has revolutionized surface science by 
cnablingtopographical measurements with atomic preci- 
sion laSSl- Despite such major advances, a difficult 
challenge that remains is the measurement of the real 
time dynamics of single molecules in their natural aque- 
ous environments. For example, the dynamics of proteins 
such as conformational changes or substrate metabolysis 
occur on force scales of 10's of piconewtons and time 
scales of 10's of milliseconds 0, 0]. This parameter 
regime is difficult to reach using current technologies. A 
promising approach is the use of micro and nanoscale 
cantilevers 111 E3, E3 - 

Atomic force microscopy relies upon detecting and in- 
terpreting the dynamics of an externally driven micron 
scale cantilever as it interacts with a sample. The stan- 
dard approaches of contact, tapping, and non-contact 
mode involve the active control and driving of the can- 
tilever probe to interact with the sample of interest. 
The ultimate sensitivity of these measurements is lim- 
ited by the underlying stochastic thermal motion of the 
cantilever. An alternative measurement approach is to 
exploit these stochastic fluctuations. This can be ac- 
complished by placing a passive (or undriven) cantilever 
in fluid and measuring the resulting stochastic dynam- 
ics. The measurement of the thermal noise spectrum is 
a commonly used calibration technique in atomic force 
microscopy [l4llT5| . 

Using a passive detection technique, successful mea- 
surement then relies upon detecting and interpreting the 
change in the stochastic dynamics of the cantilever due 
to the presence or the dynamics of the sample. For ex- 
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ample, a biomolecule could be tethered between a can- 
tilever and a surface. The stochastic dynamics of the 
cantilever would then be altered by the dynamics of the 
linking molecule and its interaction with the fluid as well 
as the proximity of the surface upon which the molecule 
is attached. The force sensitivity of this measurement is 
limited by the Brownian force on the cantilever in the 
absence of the target biomolecule, i.e. if the dynamics 
of the attached protein induce forces on the cantilever 
less than the Brownian force on the cantilever the pro- 
tein dynamics will be extremely difficult to detect. For a 
typical atomic force microscope the Brownian force on a 
single cantilever can be on the order of lOOpN (discussed 
further later) which is too l arg e for many interesting bi- 
ological measurements For example, the force 
measured by Radmacher et al. as the protein lysozyme 
metabolizes substrate is approximately 50pN 0- 

However, such a measurement technique can be sig- 
nificantly improved through the detection of correlations 
in the displacements of the two adjacent cantilevers, for 
example see Fig. QIb). The motion of one cantilever 
will cause fluid motion that will move the adjacent can- 
tilever and vice versa. The Brownian force on each can- 
tilever is uncorrelated and does not contribute to a cross- 
correlation measurement. The only contribution comes 
from the correlations due to the fluid which will be sig- 
nificantly reduced in magnitude when compared to the 
Brownian force felt by a single cantilever. For example, 
femtonewton forces have been measured from the cor- 
related motion of micron scale spheres held in closely 
spaced optical traps [TtI fl8l | . Now consider tethering a 
biomolecule between the two cantilevers in Fig. ^b). In 
this case, the motion of the two cantilevers will be cor- 
related by the fluid as well as by the dynamics of the 
linking target biomolecule. Measurements of the correla- 
tions due to the biomolecule could be used to detect the 
presence of the biomolecule and to probe its biomolecular 
dynamics in real time. 

Prior to the development of such experimental meth- 
ods the fluidic coupling between adjacent AFM can- 
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tilevers must be understood in order to build, design and 
interpret experiments where the noise from fluid induced 
correlations is small compared to the forces caused by 
the dynamics of the tethered biomolecule. The magni- 
tude of the noise induced by fluid correlations in an ar- 
ray of cantilevers is complex and depends upon many 
factors including cantilever geometry and array configu- 
ration. In this paper we explore the fluid induced cor- 
relations between two readily available atomic force mi- 
croscopes in an experimentally motivated opposing con- 
figuration. This will provide a baseline understanding 
of their stochastic dynamics in the absence of a target 
biomolecule which is essential to the success of future 
experiments. 
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FIG. 1: (a) The convention used in defining the cantilever 
geometry, not drawn to scale. Actual aspect ratios are L/h = 
98.5, w/h = 14.5, and L/w = 6.8. Specific properties of 
the cantilevers explored here are given in Table [I] (b) The 
configuration of two adjacent AFM cantilevers (not drawn to 
scale). In the absence of a tethered biomolecule between the 
cantilevers their motion is correlated due to the induced fluid 
motion. 

Although the cantilever motion is driven by molecular 
collisions the equations necessary to describe these dy- 
namics are those of classical continuum mechanics. The 
Knudscn number, Kn, represents the ratio of the mean 
free path between collisions A of a fluid molecule to a 
characteristic length L in the system. In the limit of 
Kn— > the continuum hypothesis is valid, however it 
has been found experimentally to be a good approxima- 
tion for Kn < 1 x 10 -3 |lj],[2£]. For larger values of Kn 
the continuum hypothesis fails with such consequences as 
the violation of the no-slip boundary condition at fluid 
solid interfaces. The relevant length scale L for very long 
cantilevers is the cantilever half- width w/2 from Table |I] 
and A is approximately the diameter of a water molecule, 
yielding Kn = \/L w 10~ 5 . When an array configuration 
is considered a characteristic length may also be chosen 



as the spacing between the cantilevers. For the smallest 
separation investigated here, s = 200nm, Kn w 10~ 4 . 
Since Kn is small in both situations we will assume that 
both the continuum and no-slip hypotheses are valid. 

To describe the cantilever dynamics the equations of 
elasticity l21l must be coupled with the equations of fluid 
dynamics 22] . The fluid coupling of an array of can- 
tilevers is quite complex, and as a result it is useful to 
discuss the fluid equations in more detail. The Navier- 
Stokes equations for incompressible fluids are, 

/3 ^ + R u.Vu = -Vp + V 2 m, (1) 

V«tt = 0, (2) 

where u — (u, v, w) is the fluid velocity and p is the pres- 
sure. The equations are nondimcnsionalizcd using the 
cantilever half- width w/2, the maximum cantilever oscil- 
lation velocity U, and the inverse cantilever oscillation 
frequency w _1 as the characteristic length, velocity, and 
time scales, respectively. On all material boundaries we 
impose the no-slip fluid boundary condition. The fre- 
quency parameter (3 = w 2 uj/Av is a frequency based 
Reynolds number and represents the ratio of local inertial 
forces to viscous forces where v is the kinematic viscosity 
of the liquid. The Reynolds number R = Uw/2v rep- 
resents the ratio of convective inertial forces to viscous 
forces. For the high frequency and low amplitude oscil- 
lations of interest here R <C 1 and /3 ~ 1 resulting in the 
unsteady Stokes equations. 

An often used approximation for the flow field gener- 
ated by an oscillating atomic force microscope is that 
of an infinite cylinder performing transverse oscilla- 
tions [2j|] . The fluid dynamics describing the flow field 
generated can be described in terms of viscous and poten- 
tial contributions j^]. The potential component reacts 
instantaneously and depends only upon the position of 
the cylinder and is out of phase with the cylinder mo- 
tion. The viscous component is characterized by the dif- 
fusion of momentum from the cylinder surface with dif- 
fusion constant v. The phase of the viscous component 
depends in a complicated manner upon the position of 
the cylinder. The interplay between the potential and 
viscous components of the flow field can lead to complex 
dynamics. 

An estimate of the length scale over which viscous ef- 
fects diffuse during a single cantilever oscillation is given 
by the Stokes length, 

where u>f is the resonant frequency for the beam in fluid. 
For microscale systems this length scale can become quite 
large, for example a cantilever oscillating at 50 kHz in 
water yields a Stoke's length on the order of 2 /*m. For 
very small cantilever separations, as of interest here, the 
two cantilevers will be immersed in each others Stoke's 
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layers and their dynamics will be determined by the in- 
teractions of the viscous and potential responses. It is 
important to note that although the fluid dynamics of 
simple oscillating objects such as cylinders, spheres, and 
ellipsoids are well known the coupled fluid dynamics of 
an array of such objects is poorly understood. 

The dynamics of a single micron scale cantilever in vac- 
uum can be described using the equipartition theorem for 
systems in thermal equilibrium |25(. This approach has 
been extended to include the damping effects of a viscous 
fluid [2^ . In this case, the cantilever was assumed to be 
long and thin and, as a result, the stochastic dynamics 
were described by coupling the classical elasticity equa- 
tions to the flow field caused by an oscillating infinite 
cylinder. A further simplification in [23| is the assump- 
tion that the noise force is frequency independent which 
is not strictly justified [l2Ll26| . 

A recently proposed thermodynamic approach uses the 
fluctuation-dissipation theorem to predict the stochastic 
cantilever dynamics through deterministic calculations of 
the coupled fluid-solid equations [12, Hg ■ Only a brief 
overview of the method is given here, see Refs. [HIH 
for a detailed discussion. A convenient way to use this 
approach is to calculate the deterministic dynamics of the 
cantilevers for the case where one cantilever is exposed 
to a small step force f(t) given by, 

J[l> \ fort>0. { > 

In our simulations the step force is applied to the lower 
cantilever in Fig. ^b). Upon removal of the force the 
lower cantilever returns to equilibrium through under- 
damped oscillations given by the position of the can- 
tilever tip X\ (t) . The fluid motion generated as the lower 
cantilever returns to equilibrium causes deflections in the 
tip of the upper cantilever given by X2 (t) . The auto and 
cross-correlations of the equilibrium fluctuations in can- 
tilever displacement are then given by, 

(*i(0)zi(*)> = ^Xj(t), (5) 

where Xj(t) is the deflection of the j th cantilever, 
is Boltzmann's constant, T is the absolute temperature, 
and () denotes an equilibrium average in the absence of 
the step force /. For clarity, a capital X{t) indicates 
the deterministic cantilever deflection and a small x(t) 
indicates the stochastic deflection. When j = 1 Eq. (JjjJ 
yields the autocorrelation, and for j — 2 it yields the 
cross correlation. The spectral properties of the correla- 
tions are given by the Fourier transform of the auto and 
cross-correlations to yield the noise spectra, Gn(w) and 
G12M. 



L w h k uo 



197/xm 29um 2.0/mi 1.32 N/m 452 x 10 3 rad/s 

TABLE I: Summary of the cantilever geometry and proper- 
ties. The geometry is given by the cantilever length L, width 
W, and height h. The cantilever spring constant is k and the 
resonant frequency in vacuum is o>o- The cantilevers are com- 
posed of silicon with density p c — 2320 kg/m 3 , and Youngs 
modulus E = 1.74 x 10 11 N/m 2 . The cantilevers are im- 
mersed in water with density pi = 997 kg/m 3 and dynamic 
viscosity r\ — 8.59 x 10 -4 kg/m-s. 

mj/m e uJf/uJo Q P R 



8.2 0.35 3.0 39.0 3.0 x 10" 4 



TABLE II: Summary of the cantilever dynamics in fluid as 
determined from numerical simulations in Ref [2(J: the fluid 
loaded mass of the cantilever m/, the equivalent cantilever 
mass m e , resonant frequency in vacuum coo, resonant fre- 
quency in fluid u)f, the quality Q, the frequency parameter 
/3, and the Reynolds number R. The values of m/, U)f, and 
Q have been determined by fitting the cantilever response to 
a simple harmonic oscillator. 



II. THE STOCHASTIC DYNAMICS OF AN 
ARRAY OF ATOMIC FORCE MICROSCOPE 
CANTILEVERS 

To explore the stochastic dynamics of two adjacent 
cantilevers we have performed time-dependent, three di- 
mensional finite element simulations of the governing 
deterministic fluid-solid equations (algorithm described 
elsewhere [27], E^)- The AFM cantilevers are micron 
scale with the simple beam geometry shown in Fig.^a). 
The physical properties of a cantilever are summarized 
in Table |U The stochastic dynamics of a single cantilever 
with this geometry in fluid has been explored both ex- 
perimentally |2{j and theoretically [T2I l2a | . The charac- 
teristic quantities to describe its motion are given here in 
Table [H] 26] . Most atomic force microscope cantilevers 
have a probe on the distal end that interacts with the 
sample of interest However, in this work we are in- 
terested in characterizing the fluid coupling between two 
stochastic cantilevers and, as a result, we explore the 
simpler case of two rectangular beams without attached 
probes. Experimental features such as a probe could be 
included if desired. 

A series of numerical experiments are performed for a 
range of cantilever separations s, given by 0.1 < s/h < 2 
where the separation is measured as the distance between 
the cantilever tips at zero deflection. The minimum can- 
tilever separation investigated here, s = 200nm, is a re- 
sult of computational constraints given by the numer- 
ical method used. The actual cantilever separation in 
a two-cantilever experiment will depend upon the size 
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FIG. 2: (a) The autocorrelation (solid line) and cross- 
correlation (dashed line) of the equilibrium fluctuations in 
cantilever displacement for s/h = 0.1. The two functions 
have been normalized to highlight the phase difference in the 
dynamics. The normalization used is 3 x 10 _3 nm 2 for the 
autocorrelation and 2.9 x 10 _4 nm 2 for the cross-correlation, 
(b) Cross-correlation of the equilibrium fluctuations in can- 
tilever displacement over the range of cantilever separations 
s/h = 0.1, 0.3, 0.5, 0.7, 1, 2. (inset) Detailed view of the largest 
magnitude cross-correlations illustrating a decreasing magni- 
tude with cantilever separation. 



and geometry of the cantilever probes [|| and the par- 
ticular biofunctionalization approach used to tether the 
biomolecule 0, |3(j ■ Considering readily available can- 
tilever tips and biofunctionalization protocols the range 
of cantilever separations ex plor ed spans what would be 
expected in experiment [l3Ll3l| . 

The numerical simulations yield X\(t) and X2(t) and 
the auto and cross-correlations of the equilibrium fluctua- 
tions in cantilever displacement are found using Eq. JSJ. 
The auto and cross-correlations are shown in Fig. [5] as 
a function of cantilever separation. The autocorrelation 
(xi(0)xi(t)) , normalized by its maximum value, is given 
by the solid line in Fig.[2Ia). The autocorrelation did not 



depend upon the presence of the adjacent cantilever for 
all cantilever separations tested and the curve shown is 
representative for all simulations. The cross-correlation 
for s/h = 0.1, normalized by its maximum value, is also 
plotted in Fig.^a) to clearly illustrate the phase relation- 
ship between the auto and cross-correlations. The cross- 
correlations in the equilibrium fluctuations (xi(0)x2(t)) 
are shown in Fig. |2{b). As expected, the magnitude of 
the cross-correlation decreases as the separation between 
the two cantilevers is increased (see inset of Fig. E[b)). 

It is useful to look closer at the deterministic flow field 
caused by an oscillating cantilever that would yield the 
autocorrelations through Eq. \5\ Figure |21 illustrates two 
different cross-sections of the flow field at 1.5t* where t* 
is the time at which the cantilever reaches its maximum 
velocity. This corresponds to a cantilever tip velocity 
of U = 2.5 mm/s which occurs at a time t* = 9.0/xs, 
where t*/t p ~ 0.22 and t p is the time required for the 
cantilever to complete its first oscillation. The magni- 
tude of the fluid velocity is less than 1% of the maximum 
value after a time of t rs 20i* . The numerical simulation 
used to generate the flow field of Fig. was for only a 
single cantilever in fluid to clearly illustrate the three- 
dimensional nature of the flow field around the tip of an 
oscillating cantilever. A closer inspection of the flow field 
reveals that finite fluid velocities extend further in the z- 
direction than in the y-direction as shown in panels (a) 
and (b) of Fig 0] respectively. In other words, the fluid 
flow over the cantilever tip is greater than the fluid flow 
around the sides of the cantilever. This suggests that an 
elastic object located to the side of the cantilever (in the 
y-direction) will exhibit less fluid coupling than an object 
placed at the same distance off of the tip of the cantilever 
(in the z-direction). 

Fourier transforms of the cross-correlations yield the 
noise spectra which are shown in Fig. |3Ja) as a function 
of cantilever separation. It is interesting to note that 
there is a particular frequency to* where the noise spec- 
tra vanishes, Gi2(w*) = 0. This knowledge could lead to 
experimental protocols to minimize the fluid correlated 
noise. To illustrate this further we plot the behavior of 
u>* in Fig. 13b). The solid line is a linear fit to the data 
illustrating a steady decrease in lo* with increasing sep- 
aration. The reduced frequencies for larger separations 
are a result of the finite time at which the viscous effects 
propagate. 

Using these results we can characterize the force sen- 
sitivity and time resolution of a correlated measurement 
technique using the cantilever array. An estimate of the 
force sensitivity can be found using the auto and cross- 
correlation functions j2|j , 

F U = k\{ Xl (0)xi(t))\& (6) 
F 12 = k\( Xl (0)x 2 (t))\yL' ( 7 ) 

In our notation, F\\ represents the magnitude of the 
stochastic Brownian force acting on a single cantilever. 
As shown by the solid line in Fig. Ola) the maximum 
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FIG. 3: The three-dimensional flow field near the tip (distal 
end) of an oscillating atomic force microscope. The flow field 
is a snap shot in time at the instant where the cantilever is 
at its maximum velocity, t* = 9.0/^s. (top) Flow field over 
the cantilever tip in the the x — z plane. The base of the 
cantilever is to the far left of the figure and is not shown, 
(bottom) Flow field around the sides of the cantilever in the 
x — y plane. See Fig. 0for coordinate axis definitions. In both 
figures the largest arrow indicates a fluid velocity of U = 2.5 
mm/s. 



value of the autocorrelation occurs at time t — which 
is the root-mean-squared deflection of the stochastic mo- 
tion. Using the equipartition theorem, this is given by 
the simple expression |32j. 



<* 2 > 1/2 = 



k B T 



(8) 



which yields 0.56A The magnitude of the thermally 
driven oscillations are much smaller than the thickness 
of the cantilever, (x 2 ) 1 ^ 2 /h — 2.8 x 10~ 5 . In this 
case, the Brownian force on an individual cantilever is 
F n = (k B Tk) 1/2 = 74pN. 

In a cross-correlation measurement between two can- 
tilevers the Brownian noise felt by the two individual 
cantilevers is uncorrelated and does not contribute. This 
leaves only the correlations due to the fluidic induced 
correlations. F\2 represents the approximate magni- 
tude of the force due to these hydrodynamic correla- 





y/8 



FIG. 4: (a) The variation in the x-component of the fluid ve- 
locity u as a function of distance in the z-direction measured 
from the tip of the cantilever at various times during the can- 
tilevers oscillations, (b) The variation in the x-component of 
the fluid velocity u as a function of distance in the y-direction 
measured from the tip of the cantilever at various times dur- 
ing the cantilevers oscillations. In each plot curves are shown 
for t/t* = 1, 1.5, 2, 2.5, where t/t* = 1 and 2.5 are labeled and 
the other two are in sequence. 



tions. From Fig- Etb) for the case with the closest sepa- 
ration, s/h = 0.1, the maximum magnitude of the cross- 
correlation \(xi(0)x2(t))\ max ~ 2.9 x 10~ 4 nm 2 which 
yields a root-mean-squared displacement of 0.17A Using 
Eq. @ this yields a force sensitivity of approximately 22 
pN, an improvement of better than three fold over the 
force noise acting on a single cantilever. The variation 
in F12 with cantilever separation is shown in Fig. EI The 
force due to fluid correlations decreases with separation 
and the solid line represents a quadratic polynomial curve 
fit to the data. The reduction in the noise is quite gradual 
over the separations of interest. The noise reduction at 
the largest separation s = 2/xm is a 4-fold improvement 
over a single cantilever. 

An estimate of the time resolution possible is given by 
the time it takes the cantilever to complete an oscilla- 
tion at its resonant frequency. Using the parameters ob- 
tained from the simple harmonic oscillator curve fit for 
a single cantilever in fluid, this time scale is estimated 
to be t « 2'K/ujf — 39/iS which yields a measurement 
frequency of 25 kHz. Therefore the cantilever array in- 
vestigated here can measure forces on the order of 10's 
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FIG. 5: (a) The noise spectrum Gi2(w) for the cantilever 
array, normalized by the same factor as Gn(u) in [2(|. 
Results are shown for cantilever separations of s/h = 
0.1,0.3,0.5,0.7,1,2. The curve for s/h = 2 is labelled and 
the others are in sequential order, (b) The frequency at 
which the spectral response crosses zero to* as a function 
of separation. The solid line is a linear curve fit given by 
u'/ujo = -1-58 x lQ- 2 (s/h) + 0.406. 



of pN with kHz frequency resolution which would make 
possible the real time measurements of many interesting 
molecular interactions 0, 0] . 

As mentioned, the length over which viscous effects 
act can be quite large for microscale systems. To en- 
sure that the simulation boundaries do not significantly 
affect the results a series of numerical tests were per- 
formed using varying numerical domain sizes. The size 
of the numerical domain to be used for analysis was cho- 
sen such that the magnitude of the fluid velocity went 
gradually to zero at the walls. Specifically, the distance 
between the cantilever tip and any numerical bounding 
wall is chosen such that the velocity field at t/t* = 1 
gradually decreases to a vlaue less than 1% of the max- 
imum velocity before reaching a solid wall. Numerical 
results indicate that the distance from the cantilever tip 



0.22 1 1 1 1 1 

0.5 1 1.5 2 2.5 

S/h 



FIG. 6: Force sensitivity F12 as a function of AFM cantilever 
separation. Forces have been normalized by the Brownian 
force on a single cantilever (Fu = 74 pN). The solid line is 
a quadratic curve fit to the data and is given by .F12/-F11 = 
-7.3 x lQ~ 3 (s/h) 2 - 1.33 x W~ 2 (s/h) + 0.306. 



to the walls should be at least 24<5 S in both y and z, and 
125 s in x. The numerical domain was chosen to be larger 
than these constraints in each direction for all numeri- 
cal simulations presented here. Results from simulations 
where the numerical domain was smaller than these val- 
ues yielded cantilever dynamics with lower values of Q 
and uif in agreement with the predicted increase in dis- 
sipation for a cantilever oscillating near a wall [33l l3^ | . 



III. CONCLUSIONS 

The correlated stochastic dynamics of an array of two 
AFM cantilevers have been obtained using a thermo- 
dynamic approach with deterministic numerical calcu- 
lations. Measurements of the correlations in cantilever 
displacement yield force sensitivities that are improved 
by a factor of 3 to 4 over the range of interest for can- 
tilever separations that may be useful for single molecule 
experiments. The complex flow field around the tip of 
an oscillating cantilever suggests that the configuration 
of the cantilever array is an important determining fac- 
tor in the resulting cantilever correlations. The approach 
used here is quite general and can be extended to in- 
clude more complex geometries and array configurations 
that are motivated by experiment. Theoretical calcula- 
tions, such as these, will provide important insight nec- 
essary for the interpretation and design of future micro 
and nanoscale technologies that exploit inherent thermal 
fluctuations. 
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